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Abstract 

We investigate the dependence of the longitudinal emitter dimension rn of identical bosons, 
produced in the hadronic Z° decays, on their transverse mass my obtained from 2- 
dimensional Bose-Einstein correlations (BEC) analyses. We show that this dependence is 
well described by the expression rn = cVhAt/^/mT, deduced from the uncertainty rela- 
tions, setting At to be a constant of the order of 10 -24 sec. This equation is essentially 
identical to the one previously applied to the 1-dimensional BEC results for the emitter 
radius dependence on the boson mass itself. It is further shown that a very similar be- 
haviour exists also for the dependence of the interatomic separation in Bose condensates 
on their atomic masses when they are at the same very low temperature. 
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1 Introduction 



Bose-Einstein Correlations (BEC) of identical bosons, produced in multihadron final states of 
high energy particle interactions, have been analysed for some 40 years JIJ . Many BEC analyses 
utilised pairs of identical charged pions produced in multihadron final states where the emitter 
has often been assumed to be a sphere with a Gaussian distribution. The experimental results 
have been then subjected to a Coulomb correction to account for the repulsive force between 
the equally charged particles. The kinematic variable frequently used, and still in use today, is 
defined by 

Q = \J-(qi - Q2) 2 , 

where qi and q 2 are the four momenta of the two identical hadrons. In the limit of Q — > 
the two identical bosons are occupying the same lowest energy ground state defined in their 
centre of mass system. To observe the effect of the BEC the experimental Q distribution is 
divided by a corresponding reference sample distribution which is selected so as to be, as much 
as possible, identical to the data sample in all its features but void of Bose-Einstein statistics 
effects. The distribution of this ratio is in general described by the expression 

C 2 (Q) = 1 + Xe- Q2r2 , (1) 

where r measures the average distance between the two boson when they are predominantly in 
an s-wave. This r value, extracted from C2{Q) as Q approaches zero, is taken to represent the 
dimension of the hadron emitter. The factor A in Eq. [I], which can vary between < A < 1, is 
the strength of the effect which depends on the chaoticity of the emitter and on the purity of 
the measured data sample. In e + e~ annihilations r was found to be within the range of 0.7 to 
1.0 fm (see e.g. Ref. @), essentially independent of the centre of mass energy yfsZe- 

In recent years the BEC analyses in e + e~ annihilations have been extended in several directions. 
Among them, the search for the so called higher order BEC namely, of three or more genuine 
identical hadrons correlations; the search for deviation from an ideal spherical emitter and to 
studies aimed to determine whether the correlation dimension is a function of the hadron mass. 

In Ref. || it was first pointed out that in e + e~ — > Z° — > hadrons the measured dimension 
r values are a decreasing function of the hadron mass m (see Fig. 1). This observation was 
first deduced from the measured r values obtained from BEC analyses of identical n tt and 
K ± K ± boson-pairs. That indeed dr(m)/dm < was significantly strengthened by the recent 
emitter size measurements 0] of the AA and AA pairs. These last measurements utilised a 
method proposed in Ref. || where the mixture of S=0 and S=l spin states can be determined 
for the hyperon-pair as a function of their centre of mass (CM) energy. The onset of the Pauli 
exclusion principle, as the CM kinetic energy decreases to zero and the s-wave of the systems 
dominates, determines the r{m^) value which was found to be of the order of 0.15 fm. 

Whereas the experimental findings that r(m n ) is somewhat larger than rimx), but still equal 
within errors, may still be consistent with the string fragmentation model, although in its basic 
form it expects r{m) to increase with m ||, the much smaller value obtained for r(m\) poses a 
challenge to the model p[ . At the same time however, it was shown in Ref. that by applying 
the Heisenberg uncertainty relations, one can derive an expression for r(m) which decreases as 



. . cVhAt 

rim) = =— . (2) 

\/m 



m increases, namely: 



Taking for At the value 10 -24 sec to represent the time scale of the strong interactions sec- 
tor, independent of the hadron mass, one obtains the continuous thin line in Fig. [I] which 
follows rather well the trend of the r values measured in the LEP1 data. A fit of Eq. to 
the data yields for At the value (1.2 ± 0.3) x 10~ 24 sec. The continuous thick line in Fig. [1], 
which is almost identical to the one obtained from the uncertainty relations, was derived from 
the virial theorem assuming Local Parton Hadron Duality M using a general QCD potential pfl. 
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Figure 1: The measured emitter radius r(m) as a function of the hadron mass determined from 
BEC analyses using hadronic Z° decay events at LEP1 (taken from Ref. [[|). The continuous 
thin line is the prediction from the Heisenberg uncertainty relations setting At = 10~ 24 sec, the 
upper and lower dashed lines correspond respectively to the At values of 1.5 x 10~ 24 sec and 
0.5 x 10~ 24 sec. The continuous thick line is derived from the virial theorem assuming Local 
Parton Hadron Duality and using a general QCD potential. 



The effective range of the two-pion sources was also estimated in 2-dimensional BEC analyses, 
in heavy-ion collisions |H| and in the hadronic Z° decays ||11|| , as a function of the transverse 



mass iiit of the pion-pair, defined as 



m T = 0.5 x (ym 2 + p\ T + ym 2 + p\ T j . (3) 

Here pi y and P2,t are the transverse momentum of the two identical bosons in the longitudinal 
centre of mass system (LCMS) |I2| . As can be seen e.g. in Fig. [2|, the results of these studies 
show also a decrease of the longitudinal range r z (= ru) of the tttt system as ttlt increases. 
Moreover, the behaviour of r^(%) follows very closely the dependence of r(m) which is a 
function of the hadron mass itself. 
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These findings pose the obvious question why the values of r(m) and r^(m^) essentially coincide 
when m = mr(7r). Also of interest is the question whether both or one of the two quantities, 
m and m?, are the basic variables on which the hadron emitter dimension depends on and thus 
should play an integral part in any model describing multi-hadron production. 

In the quest to understand the interrelation between r(m) and r z (m,T) we explore in Section 2 
the possibility that the dependence of r z on rax can also be described in terms of the Heisenberg 
uncertainty relations. Next in Section 3 we turn to another phenomenon related to the Bose 
statistics namely the Bose-Einstein Condensation. Here we show that at a fixed very low tem- 
perature the dependence of the interatomic atomic separation on the mass of the condensates 
atoms is proportional to l/^m atom . Finally a summary is presented in Section 4. 

2 The longitudinal dimension dependence on TTirp 

Most of the BEC analyses were carried out under the assumption that the emitter size is a 
Gaussian sphere. The possibility that the space-time extend of the particle emission region de- 
viates from a sphere and in fact is characterised by more than one dimension has been recently 
proposed [[L3|]. In particular the Lund group developed for the BEC a model based on a quan- 
tum mechanical interpretation of the string fragmentation probability ||14|| . In this model the 



correlation length in the longitudinal string direction should be larger than the corresponding 
range in the transverse direction. 

The experimental analyses, most of which were carried out with identical charged pions, have 
utilised the longitudinal centre of mass system. This coordinate system is defined for each pair 
of identical pions as the system in which the sum of the pion-pair momenta p\ + P2, referred 
to as the 'out' axis, is perpendicular to the 'thrust' (or jet) direction defined as the z-axis. 
The momentum difference of the pion-pair Q is then resolved into the longitudinal direction 
Q z = Qu parallel to the thrust axis, Q ou t is collinear with the pair momentum sum and the 
third axis Q S id e , is perpendicular to Q z and Qout- In this system the projections of the total 
momentum of the particle-pair onto the longitudinal and side directions are equal to zero. In 
particular pi z = —pi,zi where the index 1 and 2 refer to the first and second pion, so that 
Qz = Pi,z — P2,z = 2pi )2 = 2p z . The difference in the emission time of the pions couples 
to the energy difference between the particles only in the Q ou t direction. In a 2-dimensional 
analysis one defines the transverse component of Q by the relation 

Qt = Qout + Qlide ■ 

Thus the correlation function, which is fitted to the data, is of the form 

C 2 (Q z ,Qt) = 1 + Xe-^ Q ' +r T Q ^ , (4) 

where r z , estimated from Eq. |] as Q z approaches zero, is the longitudinal geometrical radius and 
Tt is a mixture of the transverse radius and the emission time. The experimental findings, both 
in heavy ion collisions |D| and in e + e~ annihilations HJ, verified the theoretical expectations 



that rT/r z is significantly smaller than one. 



In the frame work of the azimuthally symmetric sources of pair of identical bosons |T7| another 



variable of the emission function is considered namely, the transverse mass defined by Eq. 
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Figure 2: Preliminary results of DELPHI |18[ for the dependence of the longitudinal emitter 
dimension r z on the transverse mass my in hadronic Z° decays. The data is compared with 
the expression for r z given in Eq. ^ setting At to the best fitted value of 2.1 xlO -24 sec 
(continuous line) and the expectation for 1.0 xlCT 24 sec (dashed line). 



[3| Preliminary results of DELPHI ||T8[| concerning the dependence of r^mr) on rriT, measured 
for identical charged pion pairs present in the hadronic Z° decays, is shown in Fig. [|. As 
can be seen, r*(%) decreases with rnx in a very similar way to the decrease of r(m) as the 
mass m increases. In fact the continuous line in the figure, which is drawn according to Eq. |] 
replacing r(m) by rj(my) and m by m^, describes well the r z measurements using for At the 
value 2.1 x 10~ 24 sec. That this is the case is not surprising once one realises that Ar z , the 
longitudinal distance and Ap z , the difference in the longitudinal momentum of the two hadrons 
in the LCMS, are conjugate observables which obey the uncertainty principle 

Ap z Ar z = he , (5) 

Here Ap z is measured in GeV, Ar z = r z is given in fermi units and he = 0.197 GeV fm. In 
the LCMS one has 

Ap z r z = 2fiv z r z = p z r z = he 

where /i = m/2 is the reduced mass of the two identical hadrons of mass m and the longitudinal 
velocity v z of these hadrons. Thus 

n = - . (6) 

Vz 

Simultaneously we also utilise the uncertainty relation expressed in terms of energy and time 

AEAt = h , (7) 

where the energy is given in GeV and At in seconds. In as much that the total energy E of 
the two-hadron system is determined essentially only by their mass and their kinetic energy, 
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i.e. the potential energy can be neglected, and since in the LCMS \pi )Z \ = \p2,z\, one has 



E = E V m2 + p l* + p ly + = E V m ^ + P 2 * ' ( 8 ) 

i=l i=\ 

where and m 2,T are the transverse mass of the first and second hadron. As Q z decreases 
the longitudinal momentum p z vanishes so that we can, once < p\ < m 2 Tl expand the hadron 
energy E in terms of p 2 z /m 2 T and retain only the two first terms, 



2 



E = J2 m i,T* 



' + 4- « XXt + £#-. 0) 



Next we order the identical bosons so that > m 2j T and define 

mi t — m 2 t 
om T = — ! — — > , 

while 

miT + m 2 ,T 
m T = . 

Inserting these relations into Eq. ^| one gets, after few algebraic steps, that 

E = 2m T + 2 mT f . (10) 
As m\ is larger than (Sttit) 2 , we finally get 

2 

E w 2m T + ^- . (11) 
m T 

Since 2m^ is not a function of p z it may be considered to stay fixed as Q z 0, so that one has 

2 

AEAt = —At = h . (12) 
rriT 

Combining Eqs. |6| and |T^| one obtains 

. . cy/hAt „. 
r 2 (m T ) w . (13) 

This last equation is identical to the one derived in || for the dependence of emitter dimension 
on the boson mass when r(m) and m are replaced by r x (mr) and tut- A fit of Eq. O to the 



data shown in Fig. | yields At = (2.1 ±0.4) x 10~ 24 sec so that rJ(m T ) = 0.354/ ^m T (GeV) 
fm. This value is compatible with the value of At = (1.2 ± 0.3) x 10~ 24 sec obtained in 
for r(m) when one also takes into account the relatively wide spread of the 1-dimensional tttt 
BEC analyses results for r(m) obtained by the four LEP1 experiments (see e.g. Ref. ||). In 
heavy-ion collisions of S + Pb, at an energy of 200 GeV per nucleon, the longitudinal range 
r^(m,T) was also observed to be inversely proportional to the square root of itlt |17| namely, 



rj(my) ~ 2/ ^JrriT(GeV) fm. The ratio between the proportionality factor of 2.0 and 0.354 
may well be accounted for by the difference in the extend of the heavy ion target as compared 
to that of the e + e~ annihilation leading to hadronic Z° decays. 

The dependence of the transverse dimension, rj, on the two-pion transverse mass ttit has also 
been measured in the hadronic Z° decays f!9 |. Here again the transverse range was found to 



decrease as my increases. However unlike r z which is a geometrical quantity, tt is a mixture of 
the transverse radius and the emission time so that an application of the uncertainty relations 
is not straightforward. 
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3 Interatomic separation in Bose Condensates 



When bosonic atoms are cooled down, below a critical temperature Tg, the atomic wave-packets 
overlap and the equal identity of the particles becomes significant. At this temperature, these 
atoms undergo a quantum mechanical phase transition and form a Bose-Einstein (BE) conden- 
sate, a coherent cloud of atoms all occupying the same quantum mechanical state. This phe- 



nomenon, first predicted by A. Einstein in 1924/5, is a consequence of quantum statistics [20]. 



Detailed theoretical aspects of the Bose-Einstein condensation can be found in Ref. [Ell] and its 



up to date experimental situation is described in Ref. [22]. Concise summaries, aimed in par- 
ticular to the non-expert, both of the experimental situation and the theoretical background, 
can be found in Refs. |23|,|2~4||. To form Bose condensates one cools down, below the critical 



temperature Tg, extremely dilute gases so that the formation time of molecules and clusters 
in three-body collisions is slowed down to seconds or even minutes to prevent the creation of 
more familiar transitions into liquid or even solid states. 



The existence of BE Condensation was first demonstrated in 1995 by three groups |25[ in cool- 
ing down rubidium, sodium and lithium. Typical temperatures where BE condensates occur 
are in the range of 500 nK to 2 fiK with atom densities between 10 14 and 10 15 cm -3 . The 
largest sodium condensate has about 20 million atoms whereas hydrogen condensate can reach 
even one billion atoms. 

Let us consider a dilute homogeneous ideal gas of N identical bosonic atoms of spin zero, 
confined in a volume V. These atoms occupy energy levels e, handled here as a continuous 
variable, which are distributed according to the Bose-Einstein statistics. We further set the 
ground state to be eo = 0. If iVo is the number of atoms in this ground state and N ex is the 
number of atoms in the excited states then N = Nq + N ex . For a homogeneous ideal gas of 



identical bosonic atoms it can be shown [26| that at a low temperature T one has 



f2rrmkT\ 3/2 , s 

N ex = 2.612 W— — J , (14) 

where V is the volume occupied by the atoms of mass m. Since Tg is defined as the temperature 
where almost all bosons are still in excited states, we can, to a good approximation, equate N 
with N ex . That is 

N = 2S 12V(?^\' /2 , (15) 



For T < Tg one obtains from Eqs. [14] and ^ that the number of atoms N which are in the 
condensate state is, 

rp . 3/2" 

N = N-N ex = N ' 



1 - 



T 



B 



(16) 



where T is the temperature of the atoms lying at the excited energy states above the condensate 
energy level eo = 0. The atomic density of the Bose gas at very low temperatures, T /Tg « 1 
where N ^ Nq, is then given by 

N f27rmkT\ 3/2 . . 

"= F = 2 ' 612 (— J ' (17) 
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where k is the Boltzmann constant and p, the atomic density, has the dimension of L~ 3 . From 
this follows that p" 1 / 3 is the average interatomic separation in the Bose condensate. At the 
same time the thermal de Broglie wave length is equal to 

( h 2 V' 2 



Combining Eqs. [17] and [18] one has for the state of a Bose condensate the relation 

pA 3 B « 2.612 . (19) 

Thus the average interatomic distance in a Bose condensate, cIbe, is equal to 

d BE = p- 1 / 3 w AdB/1.378 . (20) 

Next we consider two different bosonic gases, having atoms with masses m\ and m 2 , which are 
cooled down to the same very low temperature T , below the critical temperature T B of each 
of them. In this case we will produce two Bose condensates with interatomic distances 

In / n 



dBE(rm) « f^[—jsr) ! »* = 1.2. (21) 
1.378 \mikT J 

From this follows that when two condensates are at the same fixed temperature T one has 

dBE{mi) \m~ 2 



d B E(m 2 ) 




(22) 



which is also the expectation of Eq. |2| for the dimension dependence on the mass of the hadron 
produced in high energy reactions provided At is fixed. Finally it is interesting to note that in 
as much that one is justified to replace in Eq. [21|, at a very low temperature, kT by AE and 
use the uncertainty relation AE = h/ At one derives the expression for r(m) as given by Eq. ^] 
multiplied by the factor v / 27r/l-378. 



In relating the condensates with the production of hadrons in high energy reactions one should 
however keep in mind that the interatomic separation proportionality to 1/ yfm does not neces- 
sarily imply that this should also be so for hadrons produced in high energy reactions. Common 
to both systems is their bosonic nature which allow all hadrons (atoms) to occupy the same 
lowest energy state. In addition the condensates are taken to be in a thermal equilibrium state. 
Among the various models proposed for the production of hadrons in high energy e + e~ and 
Nucleon-Nucleon reactions some attempts have also been made to explore the application of a 
statistical thermal-like model |27j|. However, whether this approach will eventually prevail is 



at present questionable. Finally the condensates are taken to be in a coherent state. In the 
case of the of hadrons one is able to measure r via BEC only if the chaoticity factor A in Eq. 
p] is different from zero i.e., only if the source is not 100% coherent. However so far there is 
no evidence for a dependence of r on A apart from that introduced by the correlated errors 
between Ar and AA produced by fitting Eq. [l| to the data. 



4 Summary 

It is shown that the longitudinal range of the two-pion emitter size obtained from 2-dimensional 
BEC analyses of the hadronic Z° decays as a function of the transverse mass, is well described 
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by the expectation of the Heisenberg uncertainty relations using a constant value for At of the 
order of 10~ 24 seconds. As a consequence it is not surprising that r^(m^) has essentially the 
same behaviour as the one observed experimentally for r(m) which is a function of the hadron 
mass itself. In both cases the range, in fm, is equal to (0.2 — 0.4) / ■Jmr. This kind of dependence 
on the transverse mass is also seen in heavy-ion collisions where it was found that r 1 « 2/ ^mx- 

It is interesting to note that the interatomic separation of atoms in different Bose condensates, 
having the same fixed temperature, is proportional to 1/ y/m, where here m is the atom mass. 
This behaviour is the same as found for the range between identical hadron-pair produced 
in high energy reactions when the time scale At is fixed. This similarity can be traced back 
to the close connection between the de Broglie wave length, applied to the bosonic atoms in 
condensates, and the Heisenberg uncertainty relations, used here to connect in high energy 
reactions the range between identical hadrons to their mass. 
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